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The /3 Cep/SPB star 12 Lacertae: extended mode 
identification and complex seismic modelling 
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ABSTRACT 

Results of mode identification and seismic modelling of the /3 Cep/SBP star 12 Lac- 
ertae are presented. Using data on the multi-colour photometry and radial velocity 
variations, we determine or constrain the mode degree, £, for all pulsational frequen- 
cies. Including the effects of rotation, we show that the dominant frequency, vi, is most 
likely a pure £ = 1 mode and the low frequency, va, is a dipole retrograde mode. We 
construct a set of seismic models which fit two pulsational frequencies corresponding 
to the modes £ = 0, pi and £ = 1, gi and reproduce also the complex amplitude of 
the bolometric flux variations, /, for both frequencies simultaneously. Some of these 
seismic models reproduce also the frequency i>a, as a mode 1 = 1, gi3 or gi4, and its 
empirical values of /. Moreover, it was possible to find a model fitting the six 12 Lac 
frequencies (the first five and va), only if the rotational splitting was calculated for 
a velocity of V ro t rs 75 km/s. In the next step, we check the effects of model atmo- 
spheres, opacity data, chemical mixture and opacity enhancement. Our results show 
that the OP tables are preferred and an increase of opacities in the Z— bump spoils 
the concordance of the empirical and theoretical values of /. 

Key words: stars: early-type - stars: oscillations - stars: rotation - stars: individual: 
f 2 Lac - atomic data: opacities 



1 INTRODUCTION 

The puzzle about the origin of pulsation instability in the 
B-type main sequence stars had persisted for many years. 
The explanation of this phenomenon is only 20 years old 
(Cox et al. 1992, Kiriakidis et al. 1992, Moskalik & Dziem- 
bowski 1992, Dziembowski & Pamyatnykh 1993, Gautschy 
& Saio 1993) and has been possible after re-computation of 
the opacity data, which included a huge number of transi- 
tions in the heavy element ions (Iglesias et al. 1992, Seaton 
1993) . These new opacities revealed the metal opacity bump 
around 2 ■ 10 K, which is called the Z— or Fe— bump, be- 
cause transitions within the iron M-shell dominate. Since 
then, many intriguing results have been obtained for the 
B-type pulsators, mainly thanks to data from the multi- 
site photometric and spectroscopic campaigns organized for 
several of such ob jects as well as from the spa ce-based mis- 
sions like MOST dWalker et al.ll2003ft . CoRoT i|Baglin et all 
l200rj ) or Kepler l|Koch et alj|20ich . One of the most impor- 
tant results was the discovery of high-order g-modes in early 
B-type pulsators in which typically low order p/g modes are 
observed. These variables were termed the j3 Cep/SPB stars 
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or hybrid pulsators and a few of them were detected up to 
now, e.g., v Eridani (Handler et al. 2004, Aerts et al. 2004, 
Jerzykiewicz et al. 2005), 12 Lacertae (Handler et al. 2006; 
Desmet et al. 2009) or 7 Pegasi (Handler 2009, Handler et 
al. 2009). Since the release of the first revised opacity ta- 
bles many improvements have been added in computation 
of these atomic data (Iglesias & Rogers 1996, Seaton 1996, 
2005). Moreover, the solar chemical mixture has been rede- 
termined several times (Grevesse & Noels 1993, Grevesse & 
Sauval 1998, Asplund et al. 2005, 2009). Despite of that, still 
a common problem is to get the high-order g-modes excited 
in models with masses larger than about 7 Mq. 

The hybrid pulsators give challenges as well as possibil- 
ities. Excitation of both low order p/g modes and high order 
g-modes allow potentially to scan almost the whole interior 
of a star. On the other hand, the problem with instability of 
high-order g-modes indicates that there is still something we 
do not understand or cannot adequately take into account 
in stellar physics. Usually, the opacity data are ,, blamed" for 
the inconsistencies and an artificial increase of the opacity 
coefficient was called, e.g., for the v Eri star (Pamyatnykh, 
Handler & Dziembowski 2004). 

In this paper we reanalyse the frequency spectrum of 
the hybrid B-type pulsator 12 Lacertae. In Section 2, we give 
an overview of the star. Section 3 is devoted to mode identifi- 
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cation, which we made both within the zero-rotation approx- 
imation and including some effects of rotation for frequencies 
which can be mostly affected by it. In Section 4, we present 
results of seismic modelling which we made in two steps, 
i.e., 1) fitting frequencies and 2) fitting the nonadiabatic 
/—parameter. We study the effects of model atmospheres, 
opacity data, chemical mixture and opacity enhancement. 
Conclusions are summarised in Section 5. 



2 THE VARIABLE STAR: 12 LACERTAE 

12 Lacertae (DD Lac, HR 8640, HD 214993) is a well-known 
pulsating star of B2III spectral type. Studies of the variabil- 
ity o f the star hav e begun in the early twentieth century 
when I Adams! |l912j) discovered its radial velocity variations. 
The period of these changes amounted to P = d . 193089 as 
determined by I Young (|l915l ). 

In the following years 1 2 Lac h a d bee n extensive l y stud - 
ied by many authors, e. g. lYound (|l919h . IChristiei (|1926l) . 
However, the physical interpretation of the radial veloc- 
ity changes wa s incorrec t; they thou ght that 12 Lac is a 
binary system. IStebbinsI (J1917T) and iGuthnickl | |l919l ) re- 
ported the light variations with periods of d .193 and 
d .1936, respecti vely, w h ich ag r eed w el l with the spectro- 
scopic value of lYound (|1915T ). iFathl (|1938T ) explained a 
shape of the light curve as an interference between two 



clos e per i ods of the order of .193. In his further stud- 
ies, IFathl (119471 ) found three periods; the well-known pri- 
mary d . 19308902, the secondary d . 164850 and the third 
d . 1316. However, t hese periods have not been confirmed 
by Ide Jager] |l953j), who determined the new value of 
the primary period d . 1930883 and found the secondary 
d . 197367. In 1956, one of the first worldwide observing 
campaigns was organized called The International Lacer- 
tae wee ks ide Jagerlll963r ). The analysis of these observa- 
tions by iBarnind |l963j) led to a discovery of four different 
periods; two already known: d 19308883, d 19735 8, and 
two new ones: d . 182127, 25 d .85. Ijerzvkiewiczl (119781 ) chose 
a more homogeneous sample of yellow passband data from 
this campaign and detected six periods: Pi = d . 1930754, 
P 2 = d . 19737314, P 3 = d . 18214713, P 4 = d . 1874527, 
P 5 = d .0951112 and P 6 = d .23582180. 

Sato (1979) examined observational data of 12 Lac to 
test stability of pulsational periods. Next, secular varia- 
tions of the period were searched by Ciurla (1987), but no 
variations larger than 0.08 s per century were found. Af- 
ter the great effort made during The International Lacer- 
tae weeks more photometric data were gathered and anal- 
ysed by Jarzebowski et al. (1979) and Jerzykiewicz et al. 
(1984). Smith (1980) examined Silll A 4567 line profile 
variations. His findings were generally in agreement with 
previous results of Jerzykiewicz (1978). The spectroscopic 
analysis based on high spectral and time resolution ob- 
servations had been carried out by Mathias et al. (1994). 
Their frequency analysis revealed a presence of 6 photo- 
metric frequencies in the radial velo city curve. Subsequently 
iDziembowski fe Jerzykiewicz! (J1999) investigated the possi- 
bility of existence of the equidistant triplet, v\, 1/4,, vz (corre- 
sponding to periods Pi, P4 and P3) and performed seismic 
modelling. 

In the years 2003-2004, the two multisite campaigns 
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Figure 1. The HR diagram with a p osition of the observa tional 
error box of 12 Lac as determined bv lHandler et al.l (120061 V The 
evolutionary tracks were computed for two values of metallicity, 
Z =0.010, 0.015, the hydrogen abundance of X = 0.7, initial 
rotational velocity of V ro t = 50 km/s, no overshooting from the 
convective core and the OPAL tables. Lines labeled as n = 1, 
n = 2 and an asterisk will be discussed latter in the text. 



for 12 Lac were organized. lHandler et al.l (|2006l ) reported 
results on the photometric campaign; they confirmed five 
from six already known frequencies and added six new 
independent modes, including the low frequency typical for 
Slowly Pulsating B-type stars (SPB). The frequency values 
determined from these data are: ui—5. 179034, ^2=5.066346, 
^=5.490167, ^4=5.334357, ^=4.24062, // 6 =7.40705, 
^7=5.30912, ^8=5.2162, ^=6.7023, za =5.8341 and 
^A=0-355 29 c /d. B ased on the colour photometry 
Handler et al.l (J2006J) provided the spherical degrees £ 
for the five strongest frequency peaks and constraints on 
£ for the low-amplitude frequencies. Then, iDesmet et al.l 
(2009) presented results on the spectroscopic campaign. 
They confirmed 10 independent frequencies discovered from 
photometry. The important result was that in their spec- 
troscopy the SPB-type frequency has been also detected. 
Thus, 12 Lac is another known hybrid pulsator in which 
the /3 Cep and SP B type modes are exc ited simultaneously. 
From spectroscopy, IDesmet et al.l (|2009l ) arrived at a unique 
identification of the mode degrees, £, and azimuthal orders, 
7n, for four pulsational frequencies. 

Dziembowski & Pamyatnykh (2008), based on results 
from the 2003 world-wide campaign, made seismic modelling 
of 12 Lac and showed that the rotation rate increases sig- 
nificantly towards the star centre. They addressed also a 
problem with excitation of low frequency mode va ■ 

In Fig.[T] we show the observational error box of 12 Lac 
in the HR diagram. The star is in an advanced phase of the 
core hydrogen burning. We adopted the effective tempera- 
ture, logT c ff = 4.375 ± 0.018 and lu minosity, logL/L© = 
4.18 ±0.16 from lHandler et all ||2006J ). In the literature var- 
ious determinations of metallicity are given. Niemczura & 
Daszynska-Daszkiewicz (2005) obtained [m/H]~ —0.2 ± 0.1 
from the ultraviolet IUE spe ctra, which corresp onds to 
Z « 0.012 ± 0.003, whereas iMorel et al.l i|2006l ) derived 
Z = 0.0089 ± 0.0018 from the high resolution optical spec- 
tra. Therefore, we show evolutionary tracks for two values 
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of metallicity: Z = 0.015 and Z — 0.010. Computations 
were performed using the Warsaw-New Jersey evolutionary 
code ( Pamyatnykh et al. 1998), with the OPAL equation of 
state (Rogers & N avfonovl 120021 ) . the OPAL opacity tables 
(llglesias fc Rogers! 1996 ) and the latest heavy element mix- 
ture bv lAsplund et all |2009l ). hereafter AGSS09. For lower 
temperatures, t he opacity data were supplemented with th e 
Ferguson tables iJFerguson at al.ll2005l ; ISerenelli at al.ll2009l) . 
We assumed the hydrogen abundance of X = 0.7, initial ro- 
tational velocity of Kot = 50 km/s and no overshooting 
from the convective core. Rotational velocity was chosen of 
the same magnitude as in lDesmet et al.l (|2009l ) who derived 
Kot sin i = 36 ± 2 km/s and V TO t = 49 ± 3 km/s. It means 
that 12 Lac rotates at a speed of about 10% of its break-up 
velocity (Kof ~ \/W ~ 480 km / s )- 



3 MODE IDENTIFICATION 

To identify pulsational modes detected in 12 Lac, we 
will use the most popular observables, i.e., the ampli- 
tudes and phases of the photometric and radial velocity 
variat ions. Here, we rely on the Stromgren uvy photom- 
etry ( Handler et al.l 120061 ') and the radial velocity changes 
(jDesmet et al.ll2009l ) determined as the first moment of the 
Silll 4552. 6A line. We performed mode identification in two 
steps. Firstly, we determined the degree £ ignoring all effects 
of rotation on pulsation. This is a commonly used approach 
and in the case of 12 Lac, which rotates at a speed of about 
50 km/s, seems to be justified for the most frequencies. How- 
ever, this value of rotation can already couple pulsational 
modes and affects properties of slow modes such as the fre- 
quency u\ detected in 12 Lac. Therefore, in the second step, 
we checked whether the dominant frequency, v\, can be a ro- 
tational coupled mode and apply traditional approximation 
to identify the angular numbers of the low frequency, u\. 



3.1 The zero-rotation approximation 

If all effects of rotation on stellar pulsation are neglected 
then identification of the mode degree, £, from photometric 
amplitudes and phases are independent of the inclination an- 
gle, i, intrinsic mode amplitude, e, and azimuthal order, m. 
This is also a disadvantage of disregarding rotation because 
the full geometry of mode cannot be determined. Within the 
zero-rotation approximation, the complex amplitude of the 
light variation in the x passband in the framework of linear 
theory can be written as (e.g., Daszyhska-Daszkiewicz et al. 
2002) 



A x (i) = -1.086sY e m (i,0)b x e (Dl e + D 2 , e + D't 



(1) 



where Y™ is the spherical harmonic, b x is disc averaging 
factor, and D x e , D 2 ,i and D x e are the temperature, geo- 
metrical and pressure effects, respectively. Expressions for 
these quantities can be found in Daszyhska-Daszkiewicz et 
al. (2002). 

The radial velocity variations due to pulsations are de- 
fined as the first moment of the well isolated spectral line 
and the formulae is given in Dziembowski (1977). 

To compute the theoretical values of the photometric 
amplitudes and phases two inputs are needed. The first one 



comes from models of stellar atmospheres and these are pass- 
band flux derivatives over effective temperature and grav- 
ity as well as limb darkening, included in the terms D x t 
and D% e , respectively. Here, we rely on tw o grids of stel- 
lar atmosphere models : LTE (JKurucz! 12004) and non-LTE 
IjLanz fc Hubeny||2007l ) and consider different values of the 
atmospheric metallicity, [m/H], and the microturbulent ve- 
locity, £t. The second input is connected with pulsational 
properties and this is the so-called nonadiabatic complex 
parameter / describing a ratio of the amplitude of the radia- 
tive flux perturbation to the radial displacement at the pho- 
tosphere level (e.g. Daszyhska-Daszkiewicz, Dziembowski, 
Pamyatnykh 2003). The /-parameter is incorporated in the 
temperature term, D x t . 

A comprehensive overview of possible sources of uncer- 
tainties in the theoretical values of the photometric ampli- 
tudes and p hases of early B-type pulsators has be e n stud ied 
in details bv lDaszvhska-Daszkiewicz fc Szewczukl l|201ll ). 

The value of / for a given mode can be either ob- 
tained from the linear nonadiabatic computations of stellar 
pulsations (Cugier, Dziembowski, Pamyatnykh 1994) or de- 
termined from observations (Daszyhska-Daszkiewicz et al. 
2003, 2005). Here, we applied the two approaches. 

To identify the degree I of pulsational modes of 12 Lac 
using the theoretical /-parameter, we compared the theo- 
retical and observational values of the photometric ampli- 
tude ratios and phase differences in various pairs of pass- 
bands. To this aim, we considered models from the error 
box corresponding to masses from 10 to 12Mq and used 
the nonadiabatic pulsational code of Dziembowski (1977). 
To check how our identification of £ is robust we considered 
two values of metallicity, Z — 0. 015 and 0.010, two chem i- 
cal mixtures, AGSS09 and GN93 (iGrevesse fc Noels! lll993h) 
two s ources of the opacit y ta bles, OPAL ijlglesias fc Rogers! 
1996) and OP (Scaton 2005), and two sets of stellar model 
atmospheres, LTE and non-LTE. The hydrogen abundance 
was fixed at X = 0.7 and no overshooting from the convec- 
tive core was allowed. We considered the values of £ from 
up to 6 and the frequency range appropriate to the obser- 
vational values. 

Using the theoretical values of /, we arrived at the fol- 
lowing identification. The dominant mode v\ as well as v 2 
are certainly dipoles {£ = 1). Identification of £ for v?, and v$ 
is also unique; they are quadruple modes (£ = 2). There is 
also no doubt about identification of £ for va\ certainly this 
is a radial mode (£ = 0). Identification of £ for the remaining 
frequencies is ambiguous: v§ is most likely a dipole, but it 
may be also £ = 2 or 3; vi can be t — 1, 2, 3, 5; for v$ we 
were not able to associate any values of £; fg is either £ = 3 
or 1 and vw can be £ — 1 or 2 or 3. Finally for the SPB like 
mode, va, we got £ — 1 or 6, but because of the strong aver- 
aging effects which increase with the harmonic degree, the 
£ — 6 identification is much less probable. These identifica- 
tions are independent of the metallicity, opacity data, mix- 
ture and atmosphere models. An influence of the different 
stellar parameters (within the error box) is also negligible. 

In the second approach the mode degree £ is determined 
simultaneously with the empirical values of / by fitting the 
theoretical values of the photometric amplitudes and phases 
to the observational counterparts (Daszyhska-Daszkiewicz 
et al. 2003, 2005). To this end we considered stellar parame- 
ters from the center and edges of the error box. We checked 
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Table 1. The most probable identification of £ for the 12 Lac 
frequencies obtained from two methods described in the text. In 
the last column, we put the values of m determined by Desmet 
ct al. (2009). 

frequency Ay the most probable £ m 

[c/d] [mmag] theoret. / empir. / Desmet2009 



1 

1 

0,1,2 



i/i=5. 179034 


38.1 


1 


1 


i/ 2 =5. 066346 


16.0 


1 


1 


i/ 3 =5. 490167 


11.1 


2 


2 


i/ 4 =5. 334357 


10.0 








i/ 5 =4.24062 


3.6 


2 


2 


;/6=7.40705 


2.0 


1,2,3 


1,2 


1/7=5.30912 


2.0 


1,2,3,5 


- 


i/ 8 =5.2162 


1.3 


? 


2 


1/9=6.7023 


1.3 


3,1 


3,1 


i/io=5.8341 


1.3 


1,2,3 


1,2 


i/ A =0.35529 


5.0 


1,6 


1,4 



also the effects of adopted model atmospheres, i.e., effects 
of the microturbulent velocity, £t , metallicity, [m/H] , as well 
as the non-LTE effects. All the above parameters do not 
change identification of £ although a quality of the fit ob- 
tained with the LTE and non-LTE models can differ signif- 
icantly. Usually, with the non-LTE atmospheres we got the 
better goodness of the fit. From this approach, we found a 
unique determination of £ for the six frequencies. For the first 
five modes with the largest amplitudes results are consistent 
with the determinations obtained with the theoretical values 
of the /—parameter. In the case of some other frequencies 
we succeed in getting better constraints on (,, i.e., v§ and 
i/io can be I = 1 or 2 and £ — 3 is excluded, vs has been 
identified uniquely as a quadruple. The degree of vi is inde- 
terminable because this frequency has not been detected in 
spectroscopy. The slow mode va is now i — 1 or 4. 

In Table[T] we list frequencies of 12 Lac and the most 
probable identification of the degree £ from the two ap- 
proaches. The last column contains the values of m as deter- 
mined by Desmet et al. (2009). As we can see, identifications 
of £ from the two methods are consistent. In the case of va, 
the common identification is £ = 1 and this is the most likely 
degree. Moreover, the visibility of modes with £ = 4 and 6 
is much lower. Our identifications of £ agree with those ob- 
tained by Handler et al. (2006). 

From the second approach, besides the degree £, we get 
the empirical value of /, which has a great asteroseismic 
potential, and the intrinsic mode amplitude multiplied by 
the aspect factor, i = eY(i, 0). These quantities will be dis- 
cussed in subsequent sections. 



3.2 Including effects of rotation 

Rotation causes that the 2^+1 degeneracy is lifted what 
allows to determine the full geometry of pulsational modes, 
i.e., the spherical harmonic degree, £, and azimuthal order, 
m. It also complicates mode identification because the pho- 
tometric amplitude ratios and phase differences become de- 
pendent on the inclination angle and rotational velocity. 

In the case of low order p- and g-modes, the most im- 
portant effects of rotation is mode coupling occurring for 
close frequency modes with the same values of m and i dif- 



fering by 2 (Soufi, Goupil & Dziembowski, 1998). The two 
modes are close if the difference between their frequencies is 
of the order of the rotational frequency, i.e., \vi — Vi\ ~ ^rot- 
Then, the complex photometric amplitude due to a pulsa- 
tional mode coupled by rotation is given as a sum of the 
photometric amplitudes of modes contributing to the cou- 
pling (Daszynska-Daszkiewicz et al. 2002) 



Ar P ^ = j: k a k AxAi), 



(2) 



where A(i) is given in Eq. 1. The coefficients a k results from 
linear perturbation theory of rotating stars and describe con- 
tributions of the Ifc-modes to the coup led mode. 

According to lDesmet et al.l 1J2009I ). the surface equato- 
rial rotational velocity of 12 Lac amounts to about 50 km/s. 
Adopting the stellar radius of about 7 Rq we get a surface 
rotational frequency of about z/ ro t = 0.14 c/d. As has been 
shown by Daszyiiska-Daszkiewicz et al. (2002), already ro- 
tation of this order can couple modes. 

The dominant mode of 12 Lac, which has been identified 
as the £ = 1 mode, has the photometric amplitude much 
higher than the other ones; its value of Av is about 2.5 
larger than for v-z . Because usually, it is the radial mode 
which dominates, we checked whether v\ could be 
the £ — mode coupled with the £ — 2 mode by 
rotation. Pulsational models which satisfy the rotational 
coupling conditions for the frequency v\ can be found in a 
very wide range of stellar parameters within the 12 Lac error 
box. 

In Fig. [2] we show an example of a mode coupling be- 
tween the £ = 0, pi and I = 2,gi modes with frequencies 
before the coupling of about 5.17 and 5.15 c/d, respec- 
tively. The parameters of the model are: M = 9.43Mq, 
logTeff = 4.3411 and logL/i Q = 3.975. The value of 
V TO t sin i — 36 km/s was kept constant and the values of V ro t 
are indicated at each point. The pure modes with £ =0, 1, 
2 are marked as big dots. The observational position for the 
frequency v\ is shown as an asterisk. As we can see the best 
agreement is obtained if v\ is considered as the pure £ = 1 
mode. The same conclusion was reached for other models of 
12 Lac and other pairs of modes. Therefore, in our further 
analysis we assume that v\ is the pure £ — 1 mode. 

In the case of high-order g-modes, whose frequencies are 
of the order of the rotational frequency, incorporation of the 
Coriolis force is crucial. One way to account for these effects 
of rotation is the so-called traditional approximation (e. g. 
Lee & Saio 1987, Townsend 2003). In the case of the low 
frequency of 12 Lac, va = 0.35529 c/d, the spin parameter, 
s = 2v Iot /v, amounts to 0.8-1.0 depending on the model. 
At such value of s, the photometric amplitudes and phases 
can be changed. The effects of rotation on photometric ob- 
servables for high-order g-modes excited in the SPB star 
models were studied by Townsend (2003). The formulas for 
the light and radial velocity variations, we use here, were 
given by Daszynska-Daszkiewicz, Dziembowski & Pamyat- 
nykh (2007). The semi-analytical form of the formula for the 
complex amplitude of the light variations due to pulsations 
in a high-order g-mode is as follows: 



^r d (;) = E,7^^.», 



(3) 



where 7^ . are coefficients of the expansion of the Hough func- 
tion into the series of the Legendre functions (see Daszyhska- 
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Figure 2. The diagnostic diagrams with Stromgren passbands u and j/, showing positions for coupled I = 0,pi and I = 2, gi modes 
excited in the model with M = 9.43AfQ, logT c fj = 4.3411 and log L/Lq = 3.975. Big dots indicate positions of the pure modes with 
£ = 0, 1, 2. The numbers at small dots are the values of the rotational velocity in km/s. The values of the inclination angle result from 
the condition V TO t sini = 36 km/s. An asterisk denotes an observational position of v\ with errors. 
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Figure 3. Left panel: the values of x 2 for the va frequency as a function of the rotational velocity. Only unstable modes in the model 
with M = 9.91M0, logT e ff = 4.3546 and log L/Lq = 4.039 were included. The corresponding values of the inclination angle result 
from a condition Vi- t sini = 36 km/s and are indicated at the top X axis. Right panel: the resulting amplitudes of the radial velocity 
variations calculated for the optimal values of e. The observational range of j4(V ra< j) is marked as horizontal lines. 



Daszkiewicz, Dziembowski & Pamyatnykh, 2007) and A(i) 
is given in Eq. 1. 

In order to determine the angular numbers of ua, we 
apply the same discriminant, x 2 , as proposed by Daszyriska- 
Daszkiewicz, Dziembowski & Pamyatnykh (2008). The x 2 



discriminant includes the best value of the intrinsic mode 
amplitude, e, to fit photometric amplitudes and phases in 
all passbands simultaneously. The value of x 2 is a function 
of both the inclination angle and rotational velocity. 

In the left panel of Fig. 3, we plot the dependence of 
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X 2 on the rotational velocity for unstable modes with fre- 
quencies close to va ■ Again the value of V TO t sin i = 36 km/s 
was kept constant. We show an example for model with the 
following parameters: M — 9.91A/0, logT e g = 4.3546 and 
logL/Z/0 = 4.039. Only modes which give \ 2 < 6 were 
shown. As we can see, all dipole modes are acceptable. For 
a possible discrimination of the azimuthal order, m, in the 
right panel of Fig. 3, we plotted the amplitude of the radial 
velocity variations obtained for the best value of e as a func- 
tion of Kot- The horizontal lines mark the observed value of 
^4(Kad) with errors. We can see that only the retrogarde 
mode (1=1, m = —1) has the radial velocity amplitude 
within the observational range. This identification does not 
change if other models of 12 Lac are considered. Thus, we 
assume that va is a dipole retrograde mode. 



3.3 Radial orders of pulsational modes 

Our assignment of the radial order, n, has been 
done within the zero-rotation approximation. For 
non-axisymmetric frequencies with the identified az- 
imuthal order, m, we estimated the centroid fre- 
quency, u n £, from the following equation 



^ , n£m 



Vnl 



a 



„i mu TOt , 



(4) 



where C n / is the Ledoux constant and i/ TOt is the 
rotational frequency. 

In the HR diagram (Fig. 1), we plotted the lines of con- 
stant period P = 0.18746 d, corresponding to the radial 
mode U4, assuming that 1/4 is the fundamental (n — 1) or first 
overtone mode (n — 2). Two values of metallicity, Z — 0.01 
and 0.015, were considered. Both lines are within the error 
box, but the n = 2 line is cut near the middle, for the lower 
metallicity models. It is caused by a lack of the main se- 
quence models reproducing 1/4 with masses lower than about 
11 M for Z = 0.01. 

In order to identify the radial order, n, of 1/4, we com- 
pared the empirical and theoretical values of the nonadia- 
batic /-parameter (Daszynska-Daszkiewicz, Dziembowski & 
Pamyatnykh, 2003, 2005). We chose models which are lo- 
cated on lines of constant period and are inside the observa- 
tional error box of 12 Lac. For these models, we computed 
the theoretical values of / for the frequency ^4 and com- 
pare them with the empirical counterparts. The results are 
presented in Fig.|4l where the real part of the /-parameter, 
/r, is plotted against its imaginary part, fi. The empirical 
values are marked as a box and the theoretical ones as lines 
with symbols corresponding to the three values of metallic- 
ity: Z =0.010, 0.013, 0.015. All models were computed with 
the OP opacities and no overshooting from the convective 
core. In the upper panel of Fig. [4] we assumed that ^4 is 
the fundamental mode, while in the bottom panel - the first 
overtone. 

As we can see, the agreement is achieved if V4 is the 
radial fundamental mode. Moreover, all these models are 
unstable, while in the case of the first overtone, only models 
to the right of the thin, black line labeled as 77 = 0.0 are un- 
stable. We need quite high metallicity and large masses to 
make z/4 unstable as a first overtone. With the OPAL data, 
identification of n is not so unambiguous but a requirement 
of instability and lower values of Z make the first overtone 




Figure 4. Comparison of the empirical (boxes) and theoretical 
(lines with symbols) values of / for the radial mode, 1/4, on the 
complex plane. The upper and bottom panels correspond to the 
hypothesis of the fundamental and first overtone radial mode, re- 
spectively. The theoretical values of the /-parameter were calcu- 
lated with the OP tables for three values of metallicity, Z = 0.01, 
0.015 and 0.02. All models reproduce the frequency V4 and are 
inside the error box of 12 Lac. The thin, black line in the bot- 
tom panel indicates the instability border. Models to the right 
of this line are unstable. In the upper panel, all models are un- 
stable. The arrows indicates the increasing masses and effective 
temperatures. 



also much less likely. The effect of the overshooting pa- 
rameter, Q v, on identification of n is completely neg- 
ligible. Thus, we adopted that U4 is the radial fundamental 
mode in our further analysis. Our identification of n for the 
frequency V4 confirms that of Dziembowski & Pamyatnykh 
(2008) 

Following iDesmet et al.l i|2009l ). we assumed that the 
dipole mode V2 is axisymmetric (m — 0). A survey of pul- 
sational models with different masses, M, metallicities, Z, 
and core overshooting parameters, a ov , showed, that if 1/4 
is the radial fundamental mode, then 1/2 can be only the gi 
mode. To account for the rotational splitting of v\ and 1/2, 
which belong to the £ = 1, gi triplet, Dziembowski& Pamy- 
atnykh (2008) suggested a nonuniform rotation with an in- 
crease of its rate between the envelope and core. Therefore, 
we considered rotational splitting according to the two val- 
ues of the rotational velocity: 50 and 80 km/s. The values of 
the Ledoux constant were appropriate for models of 12 Lac. 
Identification of n for the first five modes is consistent for 
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Figure 5. The overshooting parameter, a ov , as a function of 
metallicity, Z, for seismic models of 12 Lac found from fitting 
the frequencies 1/4 and V2 (the modes £ = 0, pi and £ = 1, gi, 
respectively) for the hydrogen abundance of X = 0.7 and the 
OPAL opacities. There are plotted also lines of constant mass 
(thin, black), effective temperature (black), luminosity (grey) and 
instability borders for the radial pi mode (thick, solid) and dipolc 
gi mode (thick, dotted). Models to the right of these lines are 
unstable. The grey area indicates models inside the error box. 



these two values of Kot- The v\ dipole mode was identified as 
the prograde mode (m = 1) and it is the second component 
of the gi triplet. The V3 quadruple mode is the prograde 
mode (m = 1) and it is a c omponent of the gi quintuplet. 
For the 1/5 quadruple mode, iDesmet et al.l (|2009r ) were able 
to derive only that m ^ what for most models fit to the 
g2 quintuplet. 

As we have shown in Section 3, va is the dipole retro- 
grade mode (£ — 1, m — —1). Identification of the radial 
order n for low frequency modes is very sensitive to the 
adopted value of the rotational velocity. Assuming V / I . ot =50 
km/s we got that va is a <?i5 or g\e mode and for K- O t=80 
km/s - gi3 or 314 are allowed. 

The remaining frequencies of 12 Lac do not have un- 
ambiguous identifications of the angular numbers so their 
radial orders cannot be determined as well. 



4 COMPLEX ASTEROSEISMOLOGY 

Our seismic modeling can be divided into two parts. First, 
we computed pulsational models which fit two well identi- 
fied axisymmetric pulsational frequencies. In the case of 12 
Lac, these are v$ and V2, identified as the radial fundamental 
and dipole gi modes, respectively. Then, the theoretical val- 
ues of the complex nonadiabatic parameter / corresponding 
to these two frequencies were compared with the empiri- 
cal counterparts. We termed this approach complex seismic 
model. Here, we construct complex seismic models for two 
sets of the opacity data (OPAL and OP), two sets of the 
model atmospheres (LTE and non-LTE), and two sets of 



the ch emical mixtu r e (AG SS09 and the one derived for 12 
Lac bv lMorel et~al] (|2006J )). 



4.1 Fitting frequencies 

Fitting two, well-identified frequencies, V4 and 1*2, for an as- 
sumed hydrogen abundance, X — 0.7, allow to relate the 
mass, M, and the overshooting parameter, a ov , of seismic 
models with their metallicity, Z. These relations are pre- 
sented in Fig. [5] where the overshooting parameter is plotted 
versus metallicity. Models with constant masses are marked 
with black, thin lines. Thick, solid line indicates instability 
border for the radial fundamental mode, while thick, dotted 
line - for the dipole gi mode. The instability borders of these 
two modes almost overlap and cover a very narrow range of 
metallicity, i.e., Z G (0.0105-0.0115). Models to the right of 
these lines are unstable. We marked also the lines of constant 
effective temperature and luminosity. The showed values of 
log T c ff and log L/Lq correspond to the center and edges of 
the observational error box. 

As we can see, there are a lot of seismic models of 12 
Lac fitting two frequencies, V4, and vi- An accuracy of this 
fitting amounts to 10~ 6 c/d. In general, models with a given 
value of the core overshooting have smaller masses, if they 
have larger metallicity. On the other hand, models with con- 
stant mass have increasing metallicity with decreasing over- 
shooting. These results are similar to those obtained for 7 
Peg (Walczak & Daszynska-Daszkiewicz 2010, Walczak et al. 
2013). The number of these seismic models can be limited 
if we take into account only unstable models and inside the 
observational error box (grey area). This slightly narrow the 
metallicity parameter, Z, which should be from about 0.010 
up to 0.024. Unfortunately, we are not able to constrain the 
overshooting parameter, a ov . There are seismic models with 
no overshooting from the convective core as well as with 
Qov ~ 0.5. 

In the next step, we will check whether these models re- 
produce other well identified frequencies, i.e., v^ — 5.490167 
c/d - a quadruple prograde mode (£ = 2, m = 1), v*, — 
4.24062 c/d - a quadruple prograde mode {1 = 2, m > 0), 
and va = 0.35529 c/d - the dipole retrograde mode (£ — 
1, m = — 1). The centroid values of these frequencies were 
calculated according to Eq. (12) for two values of the rota- 
tional velocity, 50 and 80 km/s. Lines of models fitting addi- 
tionally one of these three frequencies are applied in Fig. [6] 
The left and right panels include models computed with the 
rotational splitting for V I0 t = 50 and 80 km/s, respectively. 
The accuracy of the fitting of these three frequencies ranges 
from 10~ 5 to 5 • 10~ 4 c/d. 

As we can see, models computed with V I0 t — 50 km/s 
reproduce z/3, vz, and va as a gie mode, only for high val- 
ues of the overshooting parameter and quite high metal- 
licites. Models fitting va as a gis mode are nearly on 
the horizontal line and give constraints on overshooting: 
Qov G (0.15,0.31). Here, V5 can be only the quadruple sec- 
toral mode (1 = 2, m = +2). In the case of models com- 
puted with V r ot = 80 km/s, we obtained lower values of the 
overshooting parameter. The frequency va has now lower 
radial order: n — 13 for a ov € (0.05, 0.24) and n = 14 for 
Qov G (0.4, 0.5) The frequency v$ can be either m — +2 for 
a ov G (0.1, 0.38) or m = +1 for a ov « 0.5 and Z > 0.018. 

It is worth to notice that models fitting additionally 
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Figure 6. The same as in Fig. [5] but we added lines of models fitting the third frequency: 1/3 or v^ or va- The left and right panels 
correspond to models computed with the rotational splitting for Vrot = 50 and 80 km/s, respectively. The grey area indicates the 
observational error box. 
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Figure 7. The same as in the right panel of Fig. [6] but we marked models fitting the empirical values of the nonadiabatic /-parameter 
(hatched areas) of both the radial fundamental mode and dipole gi mode (labeled as f(y^ + v-i). In the left panel, we used the OPAL 
data and in the right panel the OP opacities. In both cases we assumed the AGSS09 mixture and the rotational splitting for V TO t = 80 
km/s. The LTE model atmospheres with the microturbulent velocity of ft = 2 km/s were used. 



the frequency va are usually on the nearly horizontal lines, 
i.e., they have nearly constant overshooting regardless of 
the metallicity. Thus, the high-order g-modes could be po- 
tentially good indicators of a ov but proper identification of 
n is crucial. A similar behaviour can be seen in the case 
of models fitting additionally 1/5 except for the small val- 
ues of metallicity, Z < 0.01. On the other hand the line of 
models fitting additionally the frequency vz corresponding 
to (£ = 2,m — +1, gri), has the largest slope, i.e., the largest 
dependence of a ov on Z. In general, models fitting pressure 
modes form the slopped lines on the the a ov vs. Z plane 
whereas models fitting gravity and mixed modes - nearly 
horizontal lines. 



As we can see from the right panel of Fig. [6] there are a 
few intersections of two lines which means that we have mod- 
els fitting four frequencies simultaneously: [va, V2,vz,va) or 
(v4,,V2,v 5 ,va) or (^4,^2,^3,^5). 



4.2 Fitting the /-parameters 

Each pulsational frequency is associated with the non- 
adiabatic complex parameter / which gives the ratio of the 
bolometric flux changes to the radial displacement. This 
quantity is very sensitive to properties of subphotospheric 
layers where mode driving occurs. In the case of B-type 
pulsators, / shows a strong dependence on metallicity and 
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Figure 8. The same as in Fig. [7] but we used the non-LTE model atmospheres in order to derived the empirical values of the /-parameter. 
There are also shown models which fit separately the real (vertical bars) and imaginary (horizontal bars) parts of /. 



opacity data. Thus, a comparison of the theoretical and em- 
pirical values of / may yield constraints on these quanti- 
ties (Daszyriska-Daszkiewicz, Dziembowski & Pamyatnykh, 
2005). 

In a set of seismic models fitting the frequencies v^ and 
1*2, we tried to find those which reproduce also both values of 
the /-parameter within the observational errors. Such mod- 
els exist and are marked as hatched areas in Fig.[7] on the 
Z vs. Qf ov planes and are labeled as f(yA + Vi). The left 
and right panels correspond to models computed with the 
OPAL and OP data, respectively. To derive the empirical 
values of /, we used the LTE model at mospheres wi th the 
microturbulent velocity of £t = 2 km/s (|Kuruczl l2004). 

As we can see, the hatched regions are quite small, what 
substantially reduces the number of models. The OPAL seis- 
mic models are slightly outside the observational error box 
of 12 Lac. They are cooler and less luminous. They have also 
a high overshooting parameter, <a ov ~ 0.35 — 0.55. The OP 
seismic models have smaller metallicity and overshooting. 
In this case, the area /(1/4 + 1*2) is larger and entirely inside 
the error box, suggesting a preference for the OP tables. 
In the case of the high-order g-mode, va, all OPAL mod- 
els fitting this frequency reproduce also the empirical values 
of /. The line (£ — 1, gi.4) crosses the area f(v4 + V2) for 
q ov ?a 0.5 and Z « 0.014. Also almost all OP models fitting 
the frequency va reproduce its value of /; except for models 
with the core overshooting of a ov « 0.45 and metallicity of 
Z « 0.020 - 0.025. In this case the line (i = 1, gi 3 ) crosses 
the area /(1A4 + 1*2) for a ov ~ 0.2 and Z ~ 0.011. Thus with 
both opacity OP tables, we get seismic models reproduc- 
ing the three frequencies V4, 1/2, va and their /-parameters 
simultaneously, but only the OP models are located inside 
the observational error box and have ,, reasonable" values of 
the core overshooting. In the case of vz and vs, no seismic 
model reproduces their values of / in the allowed range of 
parameters, neither with the OPAL and OP data. 

Comparing the OPAL and OP seismic models, we can 
conclude that for a given Z and a ov , models calculated with 
these two opacity tables have similar masses, effective tem- 



peratures and luminosities. Instability borders for the ra- 
dial fundamental and dipole gi modes are slightly larger for 
the OP tables. The positions of lines representing models 
fitting additionally v$ or va are quite similar, the differ- 
ence is that in the case of the OP models, there is no mode 
I = 2,m = +1,(72 for ^5. The biggest change is connected 
with vz. The line of the OP models fitting the frequency 1/3 
is shifted toward much larger metallicities. 

The empirical values of the /-parameter depend slightly 
on the model atmospheres, so we decided to check this ef- 
fect. In Fig. [SI we show models fitting the empirical values 
of / for the V4 and V2 frequencies derived with the non-LTE 
models of stellar atmospheres (Lanz & Hubeny, 2007). In 
the left and right panels we used the OPAL and OP opacity 
tables, respectively. Additionally, we marked models fitting 
separately the real, fa, (vertical bars) and imaginary, fi, 
(horizontal bars) parts of the /-parameter. As we can see, 
the overshooting parameter of models fitting // decreases 
with increasing metallicity. In the case of models fitting the 
real part of /, the dependence of a ov on Z is different. First, 
for low values of Z, the overshooting parameter increases 
rapidly with metallicity and then, for larger Z, it starts to 
decrease. With the LTE atmospheres, the shape of the re- 
gions fitting the real and imaginary part of / is similar; the 
main difference is in positions of the regions. The theoreti- 
cal values of / for V4 and 1/2 fit the empirical counterparts 
derived with the non-LTE atmospheres for slightly larger 
metallicity and lower core overshooting. 



4.3 Effects of the chemical mixture and the 
opacity enhancement 

To evaluate the effect of the adopted chemical mixture, 
we performed our seismic modeling assuming the photo- 
sp heric chemical elem ent abundances of 12 Lac as derived 
bv lMorel et all (|2006t ). In Fig.|U we show the results on the 
Qov vs. Z plane. The left and right panels correspond to 
computations with the OPAL and OP data, respectively. In 
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Figure 9. The same as in Fig. [5] but for models computed for the 12 Lac mixture. We marked models reproducing the empirical values 
of the nonadiabatic /-parameter of both V4 and V2, obtained with the LTE (horizontal bars) and non-LTE (vertical bars) atmospheres. 
In the left panel we used the OPAL data and in the right panel the OP tables. 



both cases, the empirical values of / were determined adopt- 
ing both the LTE (horizontal bars) and non-LTE (vertical 
bars) model atmospheres. 

In a comparison with computations performed with 
AGSS09 and X — 0.7, now the models fitting additionally 
the /-parameters of V4 and V2 have considerably lower values 
of the core overshooting parameter, while their metallicity 
parameter was not changed much. The mean value of aw 
for models fitting the /-parameters with the 12 Lac mixture 
is approximately 0.13 lower. The lower overshooting occurs 
for both OP and OPAL opacities as well as for the LTE and 
non-LTE atmospheres. The lines of constant masses and a 
position of a grey area corresponding to the error box were 
also moved towards lower overshooting. These changes are 
caused mainly by a high er hydrogen abund ance, X — 0.764, 
which we adopted after Mor el et al.l (|2006). Abundances of 
others elements had no significant influence on the results. 
As we can see lines of models fitting additionally one of the 
three frequencies V3, Vs, va, also follow the general trend, 
i.e., a shift towards slightly lower values of a ov . In the case 
of the OPAL data, models fitting the /-parameters are still 
mostly outside the error box. Amongst seismic models com- 
puted with the OP tables, there are those that fit also the 
frequency va and its values of the /-parameter. Thus, also 
with the 12 Lac mixture, there exist seismic models fitting 
these three frequencies and corresponding values of the /- 
parameter simultaneously. 

To summarize, we constructed seismic models 
with two sets of opacity tables, two sets of model 
atmospheres, two values of microturbulent velocity 
and two mixtures of chemical elements. Now, the 
question is which case is ,,the best". To pre-answer 
this question, in Table[2l we listed the models with 
the minimum values of \ 2 f° r each case. These mod- 
els reproduce the centroid frequencies, v-2 and 1/4, and 
their empirical values of the /-parameter. The val- 
ues of x 2 were obtained from fitting the photometric 
and radial velocity amplitudes and phases when de- 



termining the empirical values of / for vi and V4 (see 
Sect. 3.1 for the description of the method). 

In the next step, we were trying to find a model fitting 
the first five frequencies, v\, V2, V4, v-i, ^5, as well as va- We 
did it by setting the value of the rotational velocity which 
gives an adequate splitting for vi, Vz, Vs, va- We obtained 
Vrot =74.6 km/s what is consistent with the value suggested 
by Dziembowski & Pamyatnykh (2008) for the inner part of 
the star. In the case of these four mode a significant contri- 
bution to the kinetic energy density comes from the g-mode 
propagating zone and their splitting is mainly determined 
by the inner layers which presumably rotate faster. This 
rotational rate agrees also with the value of V TO t obtained 
from identification of the angular numbers of the g-mode 
frequency, va, (cf. the right panel of Fig. 3). In Fig. 1101 we 
compare the observed frequency spectrum of 12 Lac with 
the theoretical one corresponding to the model reproducing 
the six frequencies. The model has the following parameters: 
M = 10.28M©, logT cft = 4.347 and log L/L = 4.00, and 
was computed for the metallicity Z — 0.0115, overshooting 
parameter a ov = 0.39, OPAL tables and the 12 Lac mixture. 
We called it MODEL-6^ and it is marked as an asterisk in 
the HR diagram (Fig. 1) . The theoretical p-modes were split 
according to the rotational rate of Kot = 50 km/s which 
corresponds to the outer layers. 

As we can see in Fig. 1101 the dipole prograde mode 
reproduces very well the dominant frequency, V\. The fre- 
quency Vq, identified as £ = 1 or 2 is quite close to the 
centroid of £ = 1, P2 and to £ = 3, po with m — +2, +3. The 
frequency vg could be £ — 3, m — —3, po or I = 2, m — 0, 
po but £ = 2 was excluded in our identification for this fre- 
quency. The value of vj is not far from £ = 2,m = 0, gi 
and for vg, identified as £ = 2, the closest counterpart is 
£ — 3, m — +1, g2. The frequency i^io, identified as £ — 1 or 
2, is very close to £ = l,m = +1, pi. 

In Table|31 we compare the theoretical values of the /- 
parameter for MODEL-6^ with the empirical counterparts 
derived with the LTE model atmospheres. The first and sec- 
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Table 2. A list of the „best" seismic models for each considered case. The basic stellar parameters are given in columns 2 — 6. The last 
two columns contains the values of x as determined from fitting the photometric and radial velocity amplitudes and phases for U2 and 
f4, respectively. 
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0.45 
0.40 
0.30 


9.3027 
9.5440 
9.8440 


4.3463 
4.3513 
4.3501 


3.9711 
3.9973 
4.0003 


2.00 
3.54 

2.77 


2.76 
4.05 
3.45 


OPAL mixl2Lac LTE 
OPAL mixl2Lac NLTE 


0.0100 
0.0125 


0.41 
0.25 


10.289 
10.922 


4.3522 
4.3549 


4.0217 
4.0469 


3.49 
4.54 


3.84 
6.58 


OP LTE it = 2 
OP NLTE it = 2 


0.0115 
0.0130 


0.32 
0.24 


10.152 
10.462 


4.3693 
4.3704 


4.0846 
4.0962 


8.95 
6.20 


5.80 
5.32 


OP mixl2Lac LTE 
OP mixl2Lac NLTE 


0.0100 
0.0115 


0.21 
0.12 


11.494 
11.899 


4.3733 
4.3749 


4.1322 
4.1464 


10.98 
7.62 


6.46 
5.11 



3 - 



I 



2 - 



1 - 



- 



1 r 



y/i 



i) 



".A* »„ 



D„ 



D, 



m 



-3 -2 -1 
• 



m = -2 -1 1 
g 2 " ■■ 



V 2 

12 3 
** ** 



2 g 



AA e i4 



§1 



0.2 



0.6 



i). 



^ 



§i 
• * 



Po 



*■ • •*■ 



P. 



m 



-l o ; 

AAA 



Pi 



• - 



Pi 



Pf 



5 6 

v [c/d] 



7 



8 



Figure 10. Comparison of the observed frequencies of 12 Lac (vertical lines) with theoretical ones (symbols) corresponding to MODEL- 
6v fitting six frequencies: u±, U2, v%, v^, v^ and Vj\. The value of the rotational splitting was computed at V TO t = 50 km/s for p-mode 
frequencies and at V TO t = 74.6 km/s for g-mode frequencies. The mode degrees t = 0, 1,2, 3 were considered. 



ond lines for each frequency correspond to the empirical and 
theoretical values, respectively. There is no v-j because it was 
not found in the radial velocity variations. For the frequen- 
cies v§, vo, and z/io, we gave the two values of / obtained for 
mode degrees, (., as listed in Table 1. The differences in / re- 
sulting from adopting different model atmospheres are small; 
the /-parameter derived with the non-LTE atmospheres is 



usually slightly lower. In the first column of Table[3] we list 
the empirical values of the intrinsic mode amplitude, e, mul- 
tiplied by the spherical harmonic, |e| = \eY™(i, 0)|, where i 
is the inclination angle. 

For the radial mode, V4, the value of |e| is equal to |e|, 
i.e., to the root mean square of the the relative change of 
the stellar radius, dr/R, over the surface. Because we have 
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Table 3. The empirical values of |e| = |£|Y^ m (i, 0), the empirical 
and theoretical values of the /-parameter, the instability, r), and 
the value of x% f° r au detected frequencies of 12 Lac. We used 
the stellar model fitting six frequencies: V\, V2, v$, v^, v§ and vj^ 
(MODEL-6^). The empirical values of / were derived with the 
LTE model atmospheres. 



mode 


\e\ 


fa 


fi 


V 


v\,l = 1 


empir 0.01431(62) 
teor 


-8.41(37) 
-7.97 


-1.24(38) 
-0.11 


0.05 


v 2 ,£ = 1 


empir 0.00660(37) 
teor 


-7.69(44) 
-7.77 


-0.20(44) 
-0.23 


0.053 


v 3 ,£ = 2 


empir 0.00429(41) 
teor 


-10.7(1.3) 
-8.62 


-0.6(1.3) 
0.23 


0.041 


v 4 ,e = o 


empir 0.00631(22) 
teor 


-8.23(33) 

-8.17 


-0.31(33) 
0.06 


0.044 


v 5 ,£ = 2 


empir 0.00187(14) 
teor 


-4.53(75) 
-6.52 


-4.96(75) 
-0.90 


0.087 



V 6 ,£ = l 
v 6 ,£ = 2 



empir 0.00029(6) 


-20.6(4.0) 


4.4(4.0) 


— 


teor 


-11.59 


3.21 


-0.148 


empir 0.00046(6) 


-23.5(3.1) 


4.2(3.1) 


— 


teor 


-11.71 


3.21 


-0.153 



V8,£ = 2 



empir 0.00163(30) 



teor 



5.1(1.9) 

-8.17 



4.7(1.9) 
-0.07 



0.052 



u 9 ,e = i 

u Q ,£ = 3 



empir 0.00070(12) 
teor 

empir 0.00300(15) 
teor 



■6.0(1.3) 
-10.53 

13.3(1.8) 
-10.78 



3.3(1.3) 
1.96 

23.2(1.8) 
1.91 



-0.056 
-0.059 



u 10 ,e = i 

u 10 ,l = 2 



empir 0.00056(10) 
teor 

empir 0.00074(17) 
teor 



-7.2(1.3) 
-9.11 

-4.4(2.2) 
-9.20 



-0.9(1.4) 

0.68 
-3.0(2.2) 

0.66 



0.02 
0.023 



"A 



._ empir 0.00053(25) 
teor 



26.3(12.5) 
18.13 



2.9(12.0) 
5.37 



-0.631 



identification of the mode degree, £, and azimuthal order, 
m, for vi, V2, vz, vt and va as well as the inclination angle, 
i ~ 48°, IjDesmet et al J 12009). we were able to derive the 
relative change of the stellar radius also for these modes. 
As a result we obtained that |e| is equal to: 1.57% for vi, 
0.57% for i/2, 0.32% for v 3 , 0.63% for v 4 , 0.25% for v B (as- 
suming that £ — 2, m = 2) and 0.06% for va- The amplitude 
of the relative radius changes caused by radial pulsations 
of 12 Lac is similar to that of other /3 Cep stars. For ex- 
ample, |e| of the radial mode for 9 Oph is equal to 0.14% 
Daszvhska-Daszkiewicz & Walczak 2009), 



for v Eri - 1.72% 
Daszvhska-Daszkiewicz & Walczak 2010) and for 7 Peg - 
0.26% (jWalczak &; Daszvriska-Daszkiewicj[2oT(il ). The last 
column contains the value of the instability parameter 77. 
Some of these modes are stable (r\ < 0). The most nega- 
tive value is for the high-order g-mode, va, and amounts to 
—0.63, what asserts the well-known problem of excitation of 
high-order g-modes in the massive main sequence pulsators. 

The last issue we would like to examine is the effect 
of the opacity enhancement. The reason is that for 12 Lac, 
as well as for some other /3 Cephei stars, there was pos- 



tulated an increase of stellar opacities near the Z-bump 
JPamyatnvkh. Handler fc Dziembowskl l2004l : iMiglio et!~aj 
120071 : iDziembowski fe Pamvatnvkrj|2008f ). This should have 
solved the problem with instability of some modes as well 
as improve fitting frequencies. We decided to test the arti- 
ficial opacity enhancement near this opacity bump. To this 
aim, we increased the opacity coefficient in the OPAL tables 
by 50% in the metal bump. The obtained seismic models fit- 
ting two frequencies, 1/4 and V2, were not changed much. The 
masses, effective temperatures and luminosities were similar. 
As one could expect, the biggest difference is related to the 
instability. With higher opacities, all our models have un- 
stable modes £ = 0, pi and £ = 1, gi, even for Z — 0.007. ft 
turned out, that models calculated with the modified OPAL 
table mimic those computed with the original OP data. It 
is connected with the internal differences between the two 
opacity tables; the OP data provide larger opacity coefficient 
near the Z-bump. 

However, with the modified opacities we did not man- 
age to find models fitting the empirical values of the /- 
parameter, neither for V4 nor for V2- The increase of the opac- 
ity coefficient changed significantly the real part / whereas 
the imaginary part of / was almost insensitive to the opacity 
increase. Although the increase of the opacity parameter can 
solve problems with instabilities, a disagreement between 
the empirical and theoretical values of / argues against this 
change. The same result was obtained for 7 Pegasi (Walczak 
& Daszynska-Daszkiewicz 2010, Walczak et al. 2013). 



5 CONCLUSIONS 

The aim of the paper was to construct complex seismic mod- 
els of the /3 Cep/SPB star 12 Lacertae, i.e., models which fit 
simultaneously pulsational frequencies and the correspond- 
ing values of the nonadiabatic parameter /. As usually, we 
began with the mode identification from the multi-colour 
photometry and radial velocity data. Unambiguous deter- 
mination of £ was possible for the six frequencies of 12 Lac; 
for the remaining ones usually two possibilities were ob- 
tained. Moreover, for the dominant mode, v\, and the high- 
order g-mode, va, the effects of rotation on identification 
of £ were checked. As a result we got that v\ is the pure 
(non-rotational coupled) dipole mode and va is the dipole 
retrograde mode. 

We found plenty of models which reproduce the two 
centroid frequencies V4 and V2, identified as £ — 0, pi and 
£ — 1, gi, respectively, and their empirical values of the /- 
parameter simultaneously. These models are inside the ob- 
servational error box only if the OP opacity tables were used. 
Models computet with the OPAL data had too small effec- 
tive temperatures and luminosities. The next argument for 
the OP tables is that with this data, there exist models 
which fit additionally the high-order g-mode frequency, va, 
identified as t = 1, m = — 1, gi3, and its values of / for 
quite low values of the core overshooting, a ov ~ 0.2. With 
the OPAL data such models exist only for a ov « 0.5. Thus, 
we managed to find models which reproduce the three fre- 
quencies and their /-parameters simultaneously. This is the 
first star for which such a good concordance was achieved. 
Moreover, we showed that if one assumes the rotational ve- 
locity of the order of 75 km/s for the rotational splitting of 
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the gravity and mixed modes, then there exist models fitting 
the four frequencies and a model which fits the six frequen- 
cies (the first five and the low frequency va)- This value 
of the rotational velocity agrees with the result obtained 
from identification of I and m for the high-order g-mode 
frequency va- It confirms also the rotational rate suggested 
by Dziembowski & Pamyatnykh (2008) for the interior of 12 
Lac. 

Similarly to the case of the hybrid pulsator 7 Pegasi, we 
got that an increase of opacities in the metal bump spoils 
an agreement between the empirical and theoretical values 
of the non-adiabatic parameter /. Thus, our conclusion is 
the same: if modifications of the opacity tables are needed 
they have to be done in a more detailed way. 
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